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We develop a variational mean-field theory of the ferromagnetic transition in compounds like 
Lanthanum-Manganite within the framework of the Double- Exchange Model supplemented by mod- 
est disorder. We obtain analytical expressions for the transition temperature, its variation with the 
valence electron-density and its decrease with disorder. We derive an expression for the conductivity 
for both the paramagnetic and the ferromangetic metallic phases, and study its dependence on the 
temperature and magnetic field. A simple relation between the resistivity in the ferromagnetic phase 
and the spontaneous magnetization is found. Our results are in a good agreement with the experi- 
mental data on transition temperatures and resistivity in the manganite compounds with relatively 
small disorder. We comment on the effects of increased disorder. 



o 
o 



CO 



OS 



i 



X 



PACS numbers: 72.15.Gd, 75.30.Kz, 75.30.Vn 



I. INTRODUCTION 



Interest has revived recently in the perovskite manganese oxides Ai_ x B x Mn03 (where A is a trivalent and B is a 
trivalent atom), which were first investigated in the 1950'sE As the doping x and the temperature T are varied, these 
manganese oxides show a rich variety of phaseaj. Particularly interesting is the doping region 0.1 < x < 0.3, where the 
compounds undergo a transition from either insulating or very high resistance metallic, paramagnetic phase at high 
temperatures to a ferromagnetic phase at low temperatures. Near the transition, the resistivity of the compounds 
changes by orders of magnitude. The application of a strong magnetic field substantially reduces this effect, thus 
giving rise to a very large negative magnetoresistance. The physical mechanism, responsible for this behaviour, has 
been recently the subject of much discussion and controversy. It was initally suggestedEl, that the CMR in manganese 
oxides can be explained within the framework of the Double-Exchange ModelQ (DE). In this model it is assumed 
that the on-site direct repulsion U is the largest energy, followed in order by the Hund's rule energy J and the 
hybridization energy t between Mn-orbitals at neighboring sites. The basic conduction step is then the interchange 
of valence between neighboring Mn : [Mn +3 Mn +4 — > Mn +4 Mn +3 ]. The basic physical idea of the DE mechanism is 
that this electron conduction is largest when the initial and the final states are degenerate. The latter requirement 
, corresponds to an alignment of the spins of the manganese ions. In the opposite case, the conduction rate is suppressed 
• by a factor of tj J. As a result, a transition from a paramagnetic to a ferromagnetic state leads-lo a dramatic increase 
of the conductivity of the compound. Using a Dynamical Mean Field calculation (DMFT)Jj the double-exchange 
mechanism w.as successfully usedE2l for a quantitative description of the experimental data in LaSrMnC>3 compounds. 
A later studytil claimed that the agreement with the experiment found in Ref. [l^ was caused by an unphysical choise 
of the density of states (DOS), and was accidental. Bu-Ua calculation by Furukawatij with several different choises of 
the local DOS confirmed the results of his earlier workLL9. 

Subsequently a calculation carried out by the authors of Ref. [H| concluded that the double exchange model alone 
could not explain the experimental data for the manganese oxides. There were two objections: (i) that the double- 
exchange model gave a transition temperature an order of magnitude larger than experiments and (ii) that the often 
observed insulating-like resistivity (resistivity increasing with decreasing temperature) could not be explained by the 
double-exchange model. It was proposed in Ref. [l3|, that for the description of mangenese oxides, one should take into 
account a continuation to the metallic state of the Jahn- Teller distortion found for the insulating antiferromagnetic 
end-member (x w 0) in.-these compounds into the x-range of interest min some kind of dynamic fashion. As shown 
by a simple calculation^, objection (i) turns out to be due to an inadequate appreciation of the energetics of the 
double-exchange process. The transition temperature is related to the difference in the electronic cohesive energy of 
the ferromagnetic and paramagnetic phases, and is not given by the transition temperature of a spin model as-|i».Ref. 
|l3| . As regards the proposal of the effects of a possible Jahn- Teller distortion, substantial theoretical effortilfllla has 
failed to produce any results which can be compared to experiments for the resistivity. ■— . 

Meanwhile, there has been further progress experimentally. It was only recently pointed outtHI, that the manganese 
oxides at similar elecron densities show two qualitatively different types of behaviour: (i) a metal-insulator transition 
near T c , which in this case has relatively low values ~ 280 K and (ii) a metallic behaviour both below (a good metal) 
and above (incoherent metal with the absolute value of the resistivity near the Mott's limit) the critical temperature, 
which is comparably high (~ 380 K).t3 The difference appears to be the amount of disorder. This would tend to 
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remove the possibility that the Jahn- Teller effects, were they to occur, have much to do with the resistivity behavior. 
Instead the question to ask is why disorder so dramatically modifies the temperature dependence of the conductivity 
in the paramagnetic phase, while simultaneously reducing the transition temperature. The relation of the resitivity 
to the magnetization M(T) in the ferromagnetic phase also depends on disorder. For small disorder the temperature 
dependent part is proportional to the M(T) 2 , while for large disorder a much stronger dependence is found. 

It was suggested! that disorder effects-due to spin-disorder, lattice polarons due to the 30% difference in the volume 
of the Mn+ 3 and the Mn+ 4 ions as well as extrinsic disorder acting in concert might be responsible for the resisitivity 
in the paramangetic phase. The possibility suggested that spin-disorder alone my be sufficient turns out, as shown 
by recent numerical studies, not to be correct Jiltip Additional randomness due to subsitution disorder has been used 
in calculations to explaij&,the experimental dataED. Isoelectronic Lao.7- x R x Cao.3Mn03 shows enormous descrease of 
the critical temperaturecJ when R is Y compared to when R is Pr. Note that the ionic radius of La 3+ is 1.02 A, of 
Pr 3+ is 1.01 A, and of Y 3+ is 0.89 A. Note, that the substituion with Y besides changing the average bond angle 
introduces disorder. 

Also very interesting is the fact that not only does spin-disorder disappear for T <C T c) lattice disorder does as 
wellEaE3. This is evidenced by the remarkable variation of the Debye- Waller factor with temperature below and above 
T c . It is clear that spin and lattice disorder act in concert and quite unusual ways. Further that quenched lattice 
disorder generates extra lattice disorder which is annealed in the ferromagnetic phase. 

If indeed the difference in the properties of the CMR materials is caused by the effect of the substitutional disorder, 
then it might be possible to account for the main features of the behavior of the "paramagnetic-metallic" compounds 
using the "pure" double-exchange model. To address this questions one of the main objectives of the present paper. 
We also consider the effect of the substitutional disorder, and show that it leads to a substantial decrease of the critical 
temperature of the para- to ferro-magnetic transition, in agreement with the observed difference in T c in different 
CMR materials. In a future paper we hope to address the more subtle issues connected with cation and other disorder 
in the mixed- valent compounds. 

The paper is organized as follows. In the next section, we develop the variation mean field theory for the double- 
cchange Hamiltonian. This is a systematization of the ansatz used in ref. 0. We calculate the spin distribution 
function, and the critical temperature of the ferromagnetic transtition. In the third section, we study the effect of 
the substitution disorder on this phase transition. In Section IV we develop a semiclassical transport theory for the 
CMR compounds, and calculate the magnetic field- and temperature- dependence of the resistivity. We close with a 
summary and discussion of future directions. 



II. THE VARIATIONAL MEAN FIELD THEORY 



In the semiclassical limit of large .spin S of the manganese ions, the effective electron Hamiltonian in the double- 
exchange model can be expressed asQ 

H c s = ^2 *ocos(%/2)cjcj + ^2 \vic\a - ^sSBcos-di (1) 

(ij) » 

where the first sum includes hopping only between the nearest-neighbour manganese ions of different valencies, the 
angle 9ij is defined as the angle between the ion spins and Sj, Vi represents the effect of the substitutional disorder, 
B is the magnetic field, and hb is the Bohr magneton. The angle d is the angle between the spin Si and the magnetic 
field. It is important to note that the assumption of large U and J compared to t makes the charge carriers effectively 
spin-less. 

Neglecting the correlations in the orientations of the neighbour spins, we represent the free energy F of the system 
in terms of the single spin orientation distribution function Pq (f2). In the mean-field approximation, the distribution 
function depends only on the angle i? between the local spin and the external magnetic field B: 

Pn = ^-Po m (2) 
In the semiclassical limit, the corresponding spin entropy is then 

S spins = - / d§ sin P* (1?) log [P„ (0)] + S° (S) (3) 



where the function 5*g pins (S) does not depend on P§, and is related to our choise of the normalization of the spin 
distribution function J d(cos-d) P# = 1. This semiclassical approximation is discussed in detail in Appendix [a|. 
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FIG. 1. The distribution P#, for T/T c = 0.5 (red), 0.7 (green), 0.9 (blue), 0.99 (cyan) 
represent the "linear" approximation of Eq. (H), appropriate for a small magnetization. 



1.01 (magenta). The dotted line 



The calculation of the energy for a given spin distribution is more complicated. When the transfer integral between 
the neighbouring sites i and j is equal to a constant value i, and the effects of the substitution disorder can be ignored, 
the electron energy is given by 



Et 



[t\ = f de po (t; e) e (4) 

J — oo 



where p is the electron density of states (DOS) corresponding to the Hamiltonian ([!]) with no diagonal disorder 
(vi — 0) and constant transfer intergral tjj = t. To account for the effects of the substitution disorder Vi, we introduce 
an effective averaged DOS defined as 

p(t,e) = (po (t,e - v)) v (5) 

where the average is performed over the distribution of the v^s. 

To obtain the total energy, in the mean field approximation we average Et [t] over the distribution of the transfer 
integrals 

E = J dt P t {t)E t [t] (6) 

The transfer integral Uj can then be expressed in terms of the polar angles and which define the 

orientations of the corresponding spins, since they uniquely define the relative angle 9. Therefore the integration over 
t in Eq. (q) can be converted to the integration over the polar angles. Using the procedure discussed in detail in 
Appendix H, we derive the effective free energy functional, and by a direct minimization obtain the following integral 
equation for the spin distribution function: 



P tf (0) = exp 



, , #i rdfa r . r» e-n e\ \ c hbSb 

-2/ — / — / d^Bm^W / de—pltocosl- ;e - f + — cos# 



(7) 



Here the parameter ^ = £(T,B) accounts for the proper normalization of the distribution function P§. The last 
term accounts for the energy pbSB cos d of the spin, tilted at the angle ■& with respect to the direction of the 
external magnetic field. Finally, the first term in the exponential of Eq. ([?]) represent the energy of the electron gas, 
which depends on the spin distribution via the effective "local" bandwidth W ~ cos 9/2, determined by the relative 
orientation of the near spins. Note, that this term depends nontrivially on *d via the relative angle 9 — 9 [<f)\, <j>, 
This nonlinear integral equation allows a straightforward numerical solution by iterations. In Fig. |l| we plot the 
distribution for different values of the scaled temperature r = T/to and dimensionless magnetic field b = B /to- 

When the magnetization of the system is small, and the spin distribution is close to uniform (e.g. when the system 
is in the paramagnetic phase in a small external field), then the distribution function 
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hole concentration, x hole concentration, x 

FIG. 2. The dependence of the critical temperature on the electron concentration. The blue and green lines correspond to 
respectively rectangular and Gaussian DOS. The red line represents the x(l — x) dependence obtained in Ref. ^[ The circles 
are the experimantal data of Ref. The left panel shows T c in units of bandwidth W, while the absolute units for the right 
panel were calculated, assuming W = 1.8eV. 



Pa (0) = ~ + <5p* W , Spa « 1 (8) 

Expanding the exponential in the right hand side of Eq. ^ in 5p$ , and keeping the terms up to the first order in 
6p0, yields 

Spt = {T,B)cosd + 0(M 2 ) , (9) 
where M (T, B) = x (T) B + O (B 2 ) is the magnetization of the system. The susceptibility x is then given by 



where the critical temperature T c is given by 

T c = d-d sin $ cos $ / ek (p, — e) p I to cos — ; e 

^0 J-oc V 2 



(11) 



It is worthwhile noting the relation of this theory to an earlier mean-field variational description of the manganese 
oxides of Ref. ^[ There, in addition to the mean-field approximation, a specific functional form of the probability 
distribution of the angle between different spins was assumed, with the system magnetization being the variational 
parameter. This should be contrasted to the method of the present paper, when the functional form of the single-spin 
distribution function is derived variationally. The general dependence of the distribution derived here turns out to be 
quite similar to the one assumed earlier. But the results of the variational procedure developed in this section should 
be more accurate besides being on firmer ground. Another advantage of the present method is that it can be used 
for the description of the effects of the substitution disorder - something, which is hard to characterize within the 
framework of Ref. ||. 

As follows from Eq. (^Tj) , the critical temperature explicitly depends on the density of states, and the resulting value 
is in fact sensitive to the actual shape of DOS. However, this has only a marginal effect on the dependence of T c on 
the concentration x. To illustrate this behaviour, in Fig. ^ we plot the critical temperature as a function of the charge 
carrier concentration x for a rectangular (blue curve) and Gassian (red curva) densities of states. For comparison, we 
also plot the x(l — x) dependence (black line), obtained in an earlier worka, and the experimental data of Ref. M. 
The model densities of states are plotted in Fig. || and compared to the DOS pt, corresonding to the hamiltonian ([j]) 
with constant transfer integral, and no diagonal disorder. The effective bandwidth of the model desities of states is 
chosen to accurately reproduce the second moment (e 2 ). Note, how accurately the Gaussian density of states fits the 
profile of pt ■ 



4 




FIG. 3. The model rectangular (green) and Gaussian (blue) densities of states, and the actual DOS p t (e). For each model 
density of states, the bandwidth is chosen such that the second moment (e 2 ) is exact. Note, how accurately the Gaussian DOS 
fits the profile of pt- 



As follows from Fig.JJ, a reasonable choise of the bandwidth W = 1.8eV consistent with the calculations in the local 
density approximation^, leads to a good agreement with the experimental data. As explained earliem, the transition 
temperature is determined essentially by the difference in the cohesive energy of the ferromagnet and the paramagnet 
by the entropy of the paramagnet. The larger bandwidth of the ferromagnet by about 20% is the essential aspect of 
the energetics in the double-exchange problem. 

Consider now the effect of substitutional disorder. Substitutional disorder increases the electron-bandwidth for the 
paramagnet. The removal of spin-disorder is then expected to decrease the change in the bandwidth on becoming a 
ferromagnet. This is explicitly borne out by the theory here. 

Since the critical temperature is directly related to the effective DOS, it is sensitive to the substitution disorder 
in the system. Assuming the Gaussian distribution of the disorder strength Uj with the standard deviation Vq and 
Gaussian "bare" DOSc-3 p g cx exp [— e 2 / (3t 2 )], we obtain: 



di) 



de (fi — e) p [t e s (#) , e] sin ■& cos $ 



where the effective transfer integral t c g is defined by the equation 
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t 2 f 2 r>n«2 <? T/2 



(12) 



(13) 



The dependence (^|) is shown in Fig. |]for different electron concentrations. As extra disorder makes the ferromagnetic 
phase less favourable, the critical temperature goes down with an increase of Vo- 
lt might be tempting to attribute the difference in critical temperature between the "type-I" and "type-II" com- 
pounds to the effect of the substitutional disorder. In such model, an effective disorder strength of Vq ~ 0.7W, would 
fully account for only ~ 30% difference in the critical temperatures of "type-I" and "type-II" compounds. We sus- 
pect that large enough lattice disorder, in concert with lattice disorder localises electronic states in the paramagnetic 
phase. New considerations then eneter in to the determination of the transition temperature. These will be discussed 
separately. Also missing-itom the discussion above is the effect of the formation of spin-polarons which must occur in 
the paramagnetic phaseoEl They would tend to decrease T c but the number of spins in the polarons is rather small 
and only a modest numerical effect on the transition temperature is expected. They are however quite important for 
the dynamics near the transition. 



III. RESISTIVITY WITHOUT LATTICE DISORDER: SEMICLASSICAL TREATMENT 



In the mean-field approximation developed in the previous section, each spin independently fluctuates around the 
averaged value denned by the magnetization of the system. From the point of view of the semiclassical transport 
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FIG. 4. The dependence of the critical temperature on the average "disorder" vq (in units of the bandwidth Wo of the 
"clean" system v = 0). Curves of different colour represent different electron concentrations (from top to bottom: x = 0.4, 0.3, 
0.2, 0.1). The distribution of the disorder energies v% is Gaussian. 



theory, that would correspond to effective independent "scatterers" located at each point of the lattice. However, 
in the ferromagnetic phase, when the spin fluctuations are small compared to the averaged value, the corresponding 
electron mean free path may be substantially larger that the (Mn) lattice spacing. In this limit, in order to estimate 
the resistivity of the system, we can use the standard semiclassical transport theory. 

We introduce the average transfer integral t = to (cos (9 a /3/2)), so that the the corresponding unperturbed Hamil- 
tonian is defined as 



ffo = i £ *c+ C/S (14) 
and rest of H is treated as the "perturbation" 



2 

<a/3> 



<ij> 

The standard plane-wave diagonalization of Hq yields the dispersion law 

£k = i [cos (k x a) + cos (k y a) + cos (k z a)] (16) 

which describes "holes" near k = and "electrons" near e.g. k = n/a (1, 1, 1). When the Fermi energy is located near 
the bottom or near the top of the band, one can define the effective mass for the electrons and the holes respectively, 
m* = *f£. 

The kinetic equation for the electron distribution function /k is |2^, : 

~ eE 7ST = T £ 1 (k m k ' } |2 5 (ek ~ £k,) (/k ~ /k,) (17) 

k 

where /°(ek) is the equilibrium distribution function, and (k \V\ k') is the matrix element of the "perturbation" jl^). 
The kinetic equation ( p!7| ) has the solution: 

A . / . w .* (k) ^2^a (18) 

where the relaxation time is defined by the following equation: 

1 3a 3 —j /■„,/, , £k+k' 



r (k) 2n 2 H 



St 2 J dk'(l + ^+*.)j(ek-ekO 



,3 



- / dk' sin 2 (kia) S (e k - e k ') (19) 
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Assuming a uniform dispersion et = e (|k|), this expression reduces to the standard result for the transport relaxation 
time 



^-^y = J dk' Wkk' (1 - cos6»kk') S (e k - £k' 



(20) 



where the scattering rate Wkk' = 3a 3 (2ir 2 Ti) 1 SP (1 + ek+k'/ (6i)), and #kk' is the angle between the vectors k and 
k'. 

Near the top and the bottom of the band, the integrals in ( [l9|) allow a straightforward analytical evaluation. For 
example, for the holes we obtain: 

^'w-X 1 - § <*■>*)? (21) 

As we pointed out before, the semiclassical approach developed in the present section is appropriate only when the 
charge carrier mean free path I 3> a. Using (|2l|), for the ratio of the mean free path to the lattice spacing near the 
top of the band we obtain: 

f P 3 

(22) 



a SP 7T fl - § (kaf 

As follows from Eq. ([l9|) , p2]) , in the absence of substitution disorder, £/a is always greater than (3/tt) * (P/SP). The 
ratio P/SP is a monotocially decreasing function of temperature in the ferromagnetic phase, and constant above the 
T c , where P/SP = 8. Therefore, since the mean free path due to the spin disorder is substantially larger than the 
effective lattice spacing, we expect that in the relevant concentration range x ~ 0.1 — 0.3 such a "pure DE" system 
would generally show the metallic behavior . Indeed, in a typical "type-II" compound Lao.7Sr .3Mn03, the resistivity 
does show the metallic behaviour dp/dT > both above and below the transitionEJ. 

The effect of nonzero magnetization (caused either by the transtion to the ferromagnetic phase, or by external 
magnetic field) on the conductivity is twofold: first, it suppresses the fluctuations in transfer integrals thus decreasing 
the corresponding scattering rate; second, the increase of the average transfer integral caused by the magnetization 
leads to a decrease of the effective mass m* ~ 1/t. Both these factors lead to a decrease of the resistivity p. For a 
small magnetization, 

p(M) = po(l - k (M/Af max ) 2 ) (23) 

where in the effective mass approximation the coefficient n is equal to 9/5. For weakly disordered manganites, the 
resitivity indeed follows Eq. (23). As seen in the inset in Fig. 5(b) the experimental value is about 2 and slowly varies 
with the electron density. Taking into account the band nonparabolicity leads to a weak dependence of k on the 
concentration x, but does not fully account for an increase of k. The variation of the resistivity in the whole range of 
the sample magnetization < M < M max is shown in Fig. ||. 



IV. SUBSTITUTION DISORDER: THE EFFECT ON RESISTIVITY 

As we pointed out in the section || of the present paper, the 30% difference in the critical temperatures of the 
"type-F and "type-II" compounds implies, that in the "disordered" compounds the effective scattering potential is 
of the order of the electron bandwidth. In such conditions, the localization effects can become important, and the 
semiclassical treatment of the previous section is no longer appropriate. 

It has been proposed, that in the "strongly disordered" ("type-F) compounds, the ferro- to paramagnetic phase 
transition drives the metal-insulator transition. In the paramagnetic phase, the "combined effort" of the substitution 
and spin disorder is sufficient to localize the charge carriers, while in the ferromagnetic phase, due to larger electron 
bandwidth and weaker spin disorder, the mobility edge is below the Fermi energyclo. 

One might be tempted to think that this mechanism of the colossal magnetoresistance of the "disordered" man- 
ganites reduces the problem to an Anderson-type transition as a function of disorder alone, where the spin disorder 
is a function of the magnetization. This is not correct since the magnetic entropy is essential to the transition which 
occurs at a finite temperature unlike the Anderson transition which occurs at T = 0. 

An important question however is wheather the resistivity near the transition can be expressed uniquely as a 
function of magnetization. If the phase transition (with or without "diagonal" disorder) is characterized by a divergent 
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magnetic correlation length scale, this may be possible. It should be remembered however that resistivity depends 
on fluctuations at large momentum transfer. In a solid with lattice disorder, the ferromagnetic correlation length 
does not uniquely characterize the important disorder at short length scales even though it may be coupled to the 
magnetization as appears to be the case in the manganites. It is also possible that for sufficiently strong disorder, the 
ferromagnetic transition is replaced by a cross-over and there is no divergent correlation length. These are probably 
the reasons why no clear indication of scaling behavior expected in a continuous quantum phase transition were 
found in the recent esistivity measurementsEa. These considerations, however, go beyond the mean-field type theory 
developed in the present paper. 

In any case, we believe that in order to fully understand the physics underlying the colossal magnetoresistance 
in doped mangnites, one also has to take into account the lattice disorder, which is coupled to the spin disorder, 
via their influence on the charge carriers. Indeed, thfi-jSirong coupling between the spin and lattice disorder was 
recently demonstrated in two independent experiments£iJO Clearly, the effective lattice disorder has it's own nontrivial 
temperature dependence, and, being coupled to the charge carriers, therefore obviously leads to substantial deviations 
from the standard picture of the "static" Anderson metal-insulator transition. However, at this point we defer the 
further description of this effect. 



V. CONCLUSIONS 



As shown in the first part of the present paper, the phase transition from paramagnetic to ferromagnetic phase 
in relatively pure manganese oxides can be successfully described by a variational approximation on the double- 
exchange Hamiltonian. The results obtained for the critical temperature and its evolution with doping and the 
chemical composition of the compound are consistent with the experimental data. The decrease of T c with modest 
disorder is also understood. 

The resistivity of the "type-II" manganese compounds can also be successfully described using the DE model. 
We showed, that e.g. the theoretical dependence of the resistivity on the sample magnetization is in a quantitative 
agreement with the expetimental data. Our calculations show the robustness of the results to the particular choice of 
the electron density of states, as should be obvious since the transition temperature depends on the difference of the 
cohesive energy of the paramagnetic and the ferromagnetic phases. 

The principal problems of the manganites left unanswered in this paper concern the properties of the "type- 1" 
compounds and the remarkable effects of disorder in them in both the dynamic and static properties. These are also 
the more subtle problems. Especially interesting is the fact that extrinsic disorder appears to promote some additional 
disorder in the paramagnetic phase which is swept away togather with the spin-disorder in the ferromagnetic phase. 
We hope to provide an answer to these questions separately. 



APPENDIX A: THE SPIN ENTROPY IN THE SEMICLASSICAL MEAN-FIELD APPROXIMATION 

In the mean field approximation, when the spin density matrix of the whole system p s (^si 1 ^ is represented as 
product of diagonal density matrices pi = p^ (sj^J of the individual spins 



Then the total spin entropy 



is represented by 



^ P ins = Tr{p s log [p s ]} 



S spins = N (^)log[p (1) (S z ) 



(Al) 
(A2) 

(A3) 



In the semiclassical approximation S 3> 1, the summation over S z can be replaced by integration. Introducing new 
variable = arccos(5 z / S), we obtain: 



S spins = NS I dcosi? pV> (S cos 0) log |pW (Scostf) 



(A4) 
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where p^ (Scost?) is normalized as follows: 



S ,1 
1= P {1) (Sz) = S dcosi? (Scosi?) 



(A5) 



We now define the spin orientation distribution function P$ ~ p^ 1 ) (Scosi?), normalized as 

y dcostf (Scostf) = 1 



(A6) 



As follows from Eqns. (|A5|), (A6), the spin orientation distribution function 



Po = \p {1) (Scostf) (A7) 

Therefore, the semiclassical spin entropy 

S spins = ~nJ d cos P# log [P#] + N log [5] (A8) 

For example, in the paramgagnetic phase, when there are no external fields, and the spind orienataion distribution 
is uniform, P# — 1/2, the semiclassical spin entropy is equal to N\og(2S), which is consistent with the exact result 
N\og(2S + 1) for S> 1. Note, that the main contribution to the semiclassical spin entropy comes actually from the 



distribution- independent term in Eq. (A8) 



The semiclassical description, however, fails for large magnetization, when the spin system is almost completely 
polarized, and the distribution function starts to change substantially on the scale of 8$ ~ 1/5. In this case, the 



original expression, Eq. (A3), should be used for the calculation of the spin entropy. 



APPENDIX B: THE VARIATIONAL FREE ENERGY FUNCTIONAL 

In the present Apendix, we calculate the variational free energy functional for the double-exchange model. Using 
Eqns. (||),@, for the electron energy we obtain: 

E e=J — J — J roismi?i J d-d 2 sm-d 2 Pi} (i^i) Pa 2) J deep\t cosl \ ;e 1 (Bl) 

while the extra spin energy 

/•27T II PTT 

E s = -B I -p- d^isini?icostfiP^(??i) (B2) 
Jo 27r Jo 

The free energy can be obtained by the subsituting these expressions and the entropy (|3|) into the standard definition 
of the free energy 

F = E e + E s — TS (B3) 

In order to find the single-spin distribution Pg, one has to minimize the effective free energy, taking into account the 
constraints of normalization. Using the standard Lagrange multiplier method, for the effective free energy functional 
we obtain: 



2tt jx /-27T J J. 

1 / #2 



F[P#;(x,\,C] = J J J ddiwu^J d0 2 Bm02P*(0i)P4(02) J de (e - A) p 1 1 cos I - ]);e 



+ T -W 1 d&ieantiPe (i?i)log[P tf (i?i)] -B / -p- d&x sintfi costfiP* 
Jo 27T Jo Jo 27T Jo 

+ C / d$ sin t? p# (0) + F [x, A, C] (B4) 



o 
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where the "constant" F represents the spin distriibution-independent part of the free energy. Here, the Lagrange 
multiplier £ accounts for the normalization of the distribution function P§ , while the Lagrange multiplier A represents 
the constraint of having a fixed concentration of mobile electrons in the system. 

It is straightforward to show by a direct calculation, that at the extremum of the functional (B4) A = p,. This has 
a clear physical meaning - the Lagrange multiplier A corresponds to the electron num ber conservation, and therefore 
shoul be equal to the electron electrochemical potential. Replacing A by p in (B4), we finally obtain the effective 
mean field free energy functional: 
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Taking the functional derivative of (B5) with respect to P?, for the distribution we obtain: 
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Note, that the exponential on the right hand side nontrivially depends on d via the angle 9 = 9 (<j>i, <?!>, 
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Temperature, T/W Q Magnetic field, /xSB/W 

FIG. 5. The resistivity and magnetization of the sample, calculated using in the variational mean-field approximation for 
the DE model, (a) the vatiation of the resistivity (red line) and magnetization (blue line) with temperature, in the absence of 
external magnetic field; (b) the resistivity as a functionref magnetization, the inset shows the variation of the coefficient k with 
concentration x, obtained from the experimental dataEH, and the theoretical prediction in the effective mass approximation 
(dashed line); (c) the resistivity as a function of temperature, for different values of the magnetic field; (d) the resistivity as a 
function of magnetic field, for different temperatures. 
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